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Abstract 

We calculate the set of 0(a s ) corrections to the double differential decay width 
dTn/(ds\ ds2) for the process B — > X s jj originating from diagrams involving the 
electromagnetic dipole operator O7. The kinematical variables si and S2 are defined as 
s i = {Pb — <7i) 2 / m b> where ph, q\, q<i are the momenta of 6-quark and two photons. While 
the (renormalized) virtual corrections are worked exactly for a certain range of s\ and 
S2, we retain in the gluon bremsstrahlung process only the leading power w.r.t. the 
(normalized) hadronic mass S3 = (pb — qi — q2) 2 / m l in the underlying triple differential 
decay width c£T 77 / \ds1ds2ds3) . The double differential decay width, based on this 
approximation is free of infrared- and collinear singularities when combining virtual- 
and bremsstrahlung corrections. The corresponding results are obtained analytically. 
When retaining all powers in S3, the sum of virtual- and bremstrahlung corrections 
contains uncanceled 1 /e singularities (which are due to collinear photon emission from 
the s-quark) and other concepts, which go beyond perturbation theory, like parton 
fragmentation functions of a quark or a gluon into a photon, are needed which is 
beyond the scope of our paper. 



1 Introduction 



Inclusive rare 5-meson decays are known to be a unique source of indirect information about 
physics at scales of several hundred GeV. In the Standard Model (SM) all these processes 
proceed through loop diagrams and thus are relatively suppressed. In the extensions of the 
SM the contributions stemming from the diagrams with "new" particles in the loops can be 
comparable or even larger than the contribution from the SM. Thus getting experimental 
information on rare decays puts strong constraints on the extensions of the SM or can even 
lead to a disagreement with the SM predictions, providing evidence for some "new physics" . 

To make a rigorous comparison between experiment and theory, precise SM calculations 
for the (differential) decay rates are mandatory. While the branching ratios for B — y X s ^ [1] 
and B — y X s £ + £~ are known today even to next-to-next-to-leading logarithmic (NNLL) 
precision (for reviews, see [21[3]), other branching ratios, like the one for B — y Xg'yy discussed 
in this paper, are only known to leading logarithmic (LL) precision in the SM jlHT]. In 
contrast to B — y X s j, the current- current operator O2 has a non-vanishing matrix element 
for b — y 577 at order a° precision, leading to an interesting interference pattern with the 
contributions associated with the electromagnetic dipole operator O7 already at LL precision. 
As a consequence, potential new physics should be clearly visible not only in the total 
branching ratio, but also in the differential distributions. 

As the process B — y X s ^ is expected to be measured at the planned Super i?-factories in 
Japan and Italy, it is necessary to calculate the differential distributions to NLL precision in 
the SM, in order to fully exploit its potential concerning new physics. The starting point of 
our calculation is the effective Hamiltonian, obtained by integrating out the heavy particles 
in the SM, leading to 
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where we use the operator basis introduced in [8]: 

1 = {s L ^T a c L ) {c Ll ^T a b L ) , 2 
o 3 = (l£7A)E,(Wz), o 4 

o 5 = {s Lllllvlp b L ) E g (qrri p q) , e> 6 

7 = ^ m b (fi) (s L a^b R ) , O s 

The symbols T a (a = 1,8) denote the SU(3) color generators; g s and e, the strong and 
electromagnetic coupling constants. In eq. (II. 2p . m&(//) is the running 6-quark mass in the 
MS-scheme at the renormalization scale fi. As we are not interested in CP-violation effects 
in the present paper, we made use of the approximation V u bV* s V t bV t * s when writing eq. 
(II. ip . Furthermore, we also put m s = 0. 

While the Wilson coefficients Ci(fi) appearing in eq. (II. ip are known to sufficient precision 
at the low scale \i ~ since a long time (see e.g. the reviews [2l|3] and references therein), 
the matrix elements (s77|Oj|£>) and (377 g\Oi\b), which in a NLL calculation are needed 
to order g 2 s and g s , respectively, are not known yet. To calculate the (Oi, (^-interference 
contributions for the differential distributions at order a s is in many respects of similar 
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complexity as the calculation of the photon energy spectrum at order a 2 needed for the NNLL 
computation. There, the individual interference contributions, which all involve extensive 
calculations, were published in separate papers, sometimes even by two independent groups 
(see e.g. [9] and [10]). It therefore cannot be expected that the NLL results for the differential 
distributions related to B — > Xg'yy are given in a single paper. As a first step in this NLL 
enterprise, we derive in the present paper the 0(a s ) corrections to the (O7, (^-interference 
contribution to the double differential decay width dT/{dsids2) at the partonic level. The 
variables si and s 2 are defined as Sj = (pb — Qi) 2 /ml, where pb and ^ denote the four- momenta 
of the 6-quark and the two photons, respectively. 

At order a s there are contributions to dTji / \ds1ds2) with three particles (s-quark and 
two photons) and four particles (s-quark, two photons and a gluon) in the final state. These 
contributions correspond to specific cuts of the fe-quarks self-energy at order a 2 xa s , involving 
twice the operator O-j. As there are additional cuts, which contain for example only one 
photon, our observable cannot be obtained using the optical theorem, i.e., by taking the 
absorptive part of the 6-quark self-energy at three-loop. We therefore calculate the mentioned 
contributions with three and four particles in the final state individually. 

As discussed in section |2j we work out the QCD corrections to the double differential 
decay width in the kinematical range 

< s 1 < 1 ; < s 2 < 1 - Si . 

Concerning the virtual corrections, all singularities (after ultra-violet renormalization) are 
due to soft gluon exchanges and/or collinear gluon exchanges involving the s-quark. 
Concerning the bremsstrahlung corrections (restricted to the same range of Si and S2), 
there are in addition kinematical situations where collinear photons are emitted from 
the s-quark. The corresponding singularities are not canceled when combined with the 
virtual corrections, as discussed in detail in section [4j We found, however, that there are 
no singularities associated with collinear photon emission in the double differential decay 
width when only retaining the leading power w.r.t to the (normalized) hadronic mass S3 = 
(Pb — <?i — Q2) 2 /m 2 m t ne underlying triple differential distribution dYn j \ds\ds2dsj) . Our 
results of this paper are obtained within this "approximation" . When going beyond, other 
concepts which go beyond perturbation theory, like parton fragmentation functions of a 
quark or a gluon into a photon, are needed [TTj . 

Before moving to the detailed organization of our paper, we should mention that the 
inclusive double radiative process B — > X s 77 has also been explored in several extensions of 
the SM [SJI7JEE2]. Also the corresponding exclusive modes, B s — > 77 and B — > Kjj, have 
been examined before, both in the SM [3H3H2T] and in its extensions [L2|IT7t[T8|l22ti30] . We 
should add that the long-distance resonant effects were also discussed in the literature (see 
e.g. (6] and the references therein). Finally, the effects of photon emission from the spectator 
quark in the B- meson were discussed in [T3"|l3"T|. 

The remainder of this paper is organized as follows. In section [2] we work out the double 
differential distribution dT 77 /(dsids2) in leading order, i.e., without taking into account QCD 
corrections to the matrix element (s^^\0 7 \b) . We retain, however, terms up to order e 1 , with 
e being the dimensional regulator (d — 4 — 2e). Section |3] is devoted to the calculation of 
the virtual corrections of order a s to the double differential decay width. In section H] the 
corresponding gluon bremsstrahlung corrections to the double differential width are worked 
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out in the approximation where only the leading power w.r.t. the (normalized) hadronic 
mass S3 is retained at the level of the triple differential decay width dT 77 /(dsids2ds 3 ). In 
section [5] virtual- and bremsstrahlung corrections are combined and the result for the double 
differential decay width, which is free of infrared- and collinear singularities, is given in 
analytic form. In section [6] we illustrate the numerical impact of the NLL corrections and in 
section [7] we present the technical details of our calculations. The paper ends with a short 
summary in section [HJ 



2 Leading order result 



In this section we discuss the double differential decay width dTj 7 / {ds\ds2) at lowest order 
in QCD, i.e. oP s . The dimensionless variables s\ and S2 are defined everywhere in this paper 

as 

Sl = o ^ S 2 = o • I 2 - 1 ) 

At lowest order the double differential decay width is based on the diagrams shown in Fig. [TJ 
The variables s\ and S2 form a complete set of kinematically independent variables for the 

b °7 s s b °7 ,s s b p °7 s b p °7 s 




Figure 1: On the first line the diagrams defining the tree- level amplitude for b — > S77 
associated with O-? are shown. The four-momenta of the 6-quark, the s-quarks and the two 
photons are denoted by pb, p s , Qi and q2, respectively. On the second line the contribution 
to the decay width corresponding to the interference of first and second diagram is shown. 



three-body decay b — > S77. Their kinematical range is as follows: 

< sx < 1 ; < s 2 < 1 - si . 

The energies E\ and E2 in the rest-frame of the fe-quark of the two photons are related to s\ 
and S2 in a simple way: Sj = 1—2 Ei/ra^. As the energies Ei of the photons have to be away 
from zero in order to be observed, the values of S\ and S2 can be considered to be smaller 
than one. By additionally requiring s\ and S2 to be larger than zero, we exclude collinear 
photon emission from the s-quark, because 2p s qi = (p s + Qi) 2 = (Pb — I2) 2 = s 2 ml > and 
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2ps<?2 = (p s + 92) 2 = (Pb — 1i) 2 = Siffij > 0. It is also easy to implement a lower cut on 
the invariant mass squared s of the of the two photons by observing that s = (qi + q 2 ) 2 = 
1 — si — s 2 . To parametrize all the mentioned conditions in terms of one parameter c (with 
c > 0), one can proceed as suggested in [5]: 

Si > c, s 2 > c, 1 — Si — s 2 > c. (2.2) 

Applying such cuts, the relevant phase-space region in the (si, s 2 )-plane is shown by the 
shaded area in Fig. [2j Our aim in this paper is to work out the double differential decay 
width in this restricted area of the s\ and the s 2 variable also when discussing the gluon 
bremsstrahlung corrections^. In this restricted region of the phase-space, the tree-level am- 




1 s 1 

Figure 2: The relevant phase-space region for (si,s 2 ) used in this paper is shown by the 
shaded area. 



plitude is free of infrared- and collinear singularities. To exhibit the singularity structure of 
the virtual corrections discussed in the next section in a transparent way, it is useful to give 
the leading-order spectrum in d = 4 — 2e dimensions. We obtain 

drS d) _ a 2 rnl^)ml\C 7 , eff (^G^\V tb V;fQ 2 d f ji_ 
ds 1 ds 2 1024 7r 5 \m b 

with 

_ [rp + e(n +r 2 + r 3 + r 4 )] (1 - s 1 - s 2 ) 
(1 - si) 2 si (1 - s 2 ) 2 s 2 

In r we retained terms of order e 1 , while discarding terms of higher order. The individual 
pieces r , . . . , r 4 read 

r = -A8s 3 2 sl + 96s 2 2 sl - 5Qs 2 sl + 8sf + 96s^s 2 - 192s 2 s 2 + 112s 2 s 2 - 56s^ a + 
112s^si - 96s 2 si + 8s l + 8sl + 8s 2 

^^In this case, the normalized invariant mass squared s of the two photons reads s = 1 — Si — S2 + S3, where 
s 3 is the normalized hadronic mass squared. The condition 1 — Sx — s 2 >c then still eliminates two-photon 
configurations with small invariant mass. 
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3252^1 — 16s 2 Si + 16s 2 



r 2 = (A8s 3 2 sl - 9Qs 2 2 sl + 56s 2 s? - 8s? - 96sijs? + 192s|s?- 

112s 2 Si + 56s 2 si - 112s 2 si + 96s 2 si - 8si - 8s 2 - 8s 2 ) log (si) 

r 3 = r 2 (si «-» s 2 ) 

r 4 = [48s^ -96s^sf + 56s 2 s?- 8s?- 96s^ + 192s^a?- 1125 2 Si + 56a|si- 
112s 2 si + 96s 2 si - 8si - 8s 2 - 8s 2 ] log (1 - s 1 - s 2 ) 

In eq. (12 .3D the symbols m&(/z) and mj denote the mass of the fc-quark in the MS-scheme 
and in the on-shell scheme, respectively. 

In d = 4 dimensions, the leading-order spectrum (in our restricted phase-space) is ob- 
tained by simply putting e to zero, obtaining 

dT?] a 2 m 2 Mml |C 7 , e// (/i)| 2 G 2 F \V tb V t *f Q 2 d (1 - * - s 2 ) . . 

r . [z.o) 



dsi ds 2 1024 vr 5 (1 - s^s^l - s 2 ) 2 s 2 

3 Virtual corrections 

We now turn to the calculation of the virtual QCD corrections, i.e. to the contributions of 
order a s with three particles in the final state. The diagrams defining the (unrenormalized) 
virtual corrections at the amplitude level are shown on the first four lines of Fig. [3] As 
the diagrams with a self-energy insertion on the external b- and s-quark legs are taken into 
account in the renormalization process, these diagrams are not shown in Fig. [3J In order to 
get the (unrenormalized) virtual corrections dT^ ce / (dsids 2 ) of order a s to the decay width, 
we have to work out the interference of the diagrams on the first four lines in Fig. [3] with the 
leading order diagrams in Fig. [TJ One of these interference contributions is shown on the last 
line in Fig. [3j To illustrate the calculational procedure for getting the virtual corrections to 
the decay width, we describe in section 17.11 the relevant steps for the particular interference 
shown in Fig. [31 

In addition, we have to work out the counterterm contributions to the decay width. They 
can be split into two parts, according to 

+ ^r- (3-1) 



dsids 2 dsids 2 dsids 2 
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Figure 3: On the first four lines the diagrams defining the one-loop amplitude for b — > 577 
associated with O7 are shown. Diagrams with self-energy insertions on the external quark- 
legs are not shown. On the last line the contribution to the decay width corresponding to 
the interference of the first diagram on the second line with the second (tree-level) diagram 
in Fig. [I] is shown. 



Part (A) involves the LSZ factors y/ and a/ Z£ s for the b- and s-quark field, as well 
as the self-renormalization constant Z^ s of the operator O7 and Z^ renormalizing the 
factor rhb(n) present in the operator O-j. The explicit results for these Z-i actors are given 
to relevant precision in Appendix O For part (A) we get 



dY 



ct,(A) 
77 



dsidsi 



5Zg s + 5Zg s + 25Z™ + 25Z™ 



dT 



(0,d) 
77 



- dsids 2 



(3.2) 



where dT^j^ / \ds1ds2) is the leading order double differential decay width in d- dimensions, 
as given in eq. (12. 3ft . 

The counterterms defining part (B) are due to the insertion of —idrrifjbb in the internal 



Figure 4: Counterterm diagrams with a 5m b insertion, see text. 



6-quark line in the leading order diagrams as indicated in Fig. HJ where 
6m b = (Z™-l)m b . 

More precisely, Part (B) consists of the interference of the diagrams in Fig.|4]with the leading 
order diagrams in Fig. [TJ 

By adding dF 77 Te / \ds\ds 2 ) and dVf 7 / (ds\ds2) , we get the result for the renormalized 
virtual corrections to the spectrum, dT 77 ,mrt /{ds\ ds 2 ). It is useful to decompose this result 
into two pieces, 

^p(l),«irt ^j-^(l,a),virt ^p(l,6),i;irt 



ds\ ds 2 dsi ds 2 ~^ ds\ ds 2 ^ 

The infrared- and collinear singularities are completely contained in dT 77 mrt / (dsi ds 2 ) . 
Explicitly, we obtain 

dr 7 Y lmrt = ^ Cf (_2 + 41og(s 1 + s 2 )-5 \ / A \ 2e rfFj^ 
ds\ds 2 An \ e 2 e J \m b J dsids 2 

where dT 77 / (dsi ds 2 ) is understood to be taken exactly as given in eqs. f 12 . 3 j) and 02 .4p . i.e., 
by including the terms of order e 1 in r. From the explicit expression drft^/idsidsj we 
see that the singularity structure consists of a simple singular factor multiplying the corre- 
sponding tree-level decay width in d- dimensions. We stress that singularities (represented by 
1/e 2 and 1/e poles) are entirely due to soft- and/or collinear gluon exchange. The infrared 
finite piece dT 77 ' mrt / (dsi ds 2 ) can be written as 



ds 1 ds 2 1024 7r 



5 



x 




— 4r n (1 — Si — So) , U y^ 2 °i Vj 

V log-^- + — 7? l= \ ^ (3.5) 



- si) 2 si (1 - s 2 ) 2 s 2 m b 3(1 - si) 3 si (1 - s 2 ) 3 s 2 
where the individual quantites v±, . . . ,v 2 o are relegated to Appendix lAl 



4 Bremsstrahlung corrections 

We now turn to the calculation of the bremsstrahlung QCD corrections, i.e. to the contribu- 
tions of order a s with four particles in the final state. Before going into details, we mention 
that the kinematical range of the variables Si and s 2 defined in eq. (12. ip is given in this case 
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by < Si < 1 ; < s 2 < 1. Nevertheless, we consider in this paper only the range which is 
also accessible to the three-body decay b — >■ 577, i.e., < Si < 1 ; < s 2 < 1 — s 1 or, more 
precisely, by its restricted version specified in eq. (I2.2p . 

The diagrams defining the bremsstrahlung corrections at the amplitude level are shown 
in the first line of Fig. [5j The amplitude squared, needed to get the (double differential) 




/ 

/ 



Figure 5: On the first line the diagrams defining the gluon-bremsstrahlung corrections to 
b — > S77 are shown at the amplitude level. The crosses in the graphs stand for the possi- 
ble emission places of the gluon. On the second line the contribution to the decay width 
corresponding to the interference of diagram 1 with diagram 6 is illustrated. 

decay width, can be written as a sum of interferences of the different diagrams on the first 
line in Fig. |5j One such interference is shown on the second line of the same figure. The 
four particle final state is described by five independent kinematical variables. In a first 
attempt we worked out the decay width by keeping si and s 2 differential and integrating 
over the three remaining variables. Proceeding in this way, we found that the infrared- 
and collinear singularities in the bremsstrahlung spectrum do not cancel when adding the 
virtual corrections. The sum still contains 1/e-poles, but no l/e 2 -poles. While, as already 
mentioned in section [3j the only source of the singularities in the virtual corrections in our 
restricted range of Si and s 2 are due to soft gluon-emission and/or collinear emission of 
gluons from the s-quark, we found after analyzing the bremsstrahlung kinematics more 
carefully that there are situations where one of the photons can become collinear with the 
s-quark. This is the reason why there is no cancellation of singularities when combining 
virtual- and bremsstrahlung corrections. Realizing that for (formally) zero hadronic mass 
of the (s, g)-system collinear photon emission is kinematically impossible, led us to the 
idea that we should first look at the triple differential decay width dTrr/(dsids 2 ds3), where 
s 3 = {Ps +Pg) 2 / m l i s the normalized hadronic mass squared. Our conjecture was that the 
double differential decay width, based on the triple differential decay width in which only 
the leading power terms w.r.t. S3 are retained, should lead to a finite result when combined 
with the virtual corrections. 

We therefore worked out the leading power of this quantity w.r.t S3, denoting it by 
c?r 1 7 c 7 admgpower /((isi(is2(is3). The leading power, which is of order I/S3 (modified by epsilon- 
tic dimensional regulators), is supposed to be a good approximation for low values of the 
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hadronic mass. An approximation to the double-differential decay width dT 77 rems / \ds\ds2) 
due to gluon bremsstrahlung corrections is then obtained by integrating 
dr 7 °j dmgpower /(dsids2ds3) over S3, which runs in the range S3 G [0, Si-s 2 ]- The approximation 
is obviously accurate for small values of s\ ■ s 2 . As s\ ■ s 2 is at most 1/4, the approximation is 
expected not to be bad in the full region of si and S2 considered in this paper. The technical 
details of the calculation of the leading power w.r.t. S3 in the triple differential decay width 
are illustrated in section 17.21 for the interference of diagram 1 with diagram 6, as shown in 
the second line of Fig. [5J 

Indeed, we find that the infrared- and collinear singularities cancel when combining the 
approximated version of dT$' brems / \ds\ds2) with the virtual corrections dT^' mrt j \ds\ds-i). 

When going beyond this approximation other concepts, which go beyond perturbation 
theory, like parton fragmentation functions of a quark or a gluon into a photon, are needed 
[TT] . We do not enter this issue in this paper. 

The result of combined virtual- and bremsstrahlung corrections is explicitly presented in 
the next section. 



5 Final result for the decay width at order a s 

The complete order a s correction to the double differential decay width dTn/ (ds\ ds2) is ob- 
tained by adding the renormalized virtual corrections from section [3] and the bremsstrahlung 
corrections discussed in section HI Explicitly we get 



dT 



(i) 
77 



a 2 m 2 M ml |C 7 , e// (/i)| 2 G\ \V tb V t * s \ 2 Q\ 



ds\ d$2 



Os 
4lT 



C f 



1024 vr 5 
-4r (1 - si - s 2 



si ) 2 s 1 (l 



log— + / 

s 2 Y s 2 m b 



where / is decomposed as 

(1 - si - s 2 ) (/i • h ■ h ■ I: ■ h ■ ./o • ./V • /s • • /ib + /i6 + hi) 



f 



3(l-si) 

ho + hi + h2 + fi3 + fu 

3(1 - sifsi (1 - s 2 ) 3 s 2 



8i(l 



\3 
82) s 2 



The individual quantities hi ■ ■ ■ > fvr read 



fx 



-16tt 2 s3s^ - 16n z s 2 s-{ + 48vr z s^s° + 48^5^ - 48vr z s 2 s° - 48vr z s°Si + 
16vr 2 s^ + 16vr 2 s^ - 168s? - 168s^ + (1104 + 160rr 2 ) s A 2 s{ + 
(-3360 - 400rr 2 ) s\s\ + (3432 + 304vr 2 ) sls\ + (-1296 - 64vr 2 ) s 2 s\ + 
(120 - 16rr 2 ) s\ + (-3360 - 400rr 2 ) s^s? + (10416 + 11527T 2 ) sls\ + 
(-10872 - 1056rr 2 ) sjsl + (3984 + 368vr 2 ) s 2 s? + (3432 + 304vr 2 ) s^s 2 + 
(-10872 - 1056rr 2 ) s 3 2 sj + (12096 + 1088rr 2 ) sjsj + 
(-4872 - 448vr 2 ) s 2 s 2 + (216 + 16tt 2 ) s\ + (-1296 - 64rr 2 ) s\ Sl + 



.2 „5„3 



2„5„2 



r 2„5. 



(5.1) 



(5.2) 
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(3984 + 368tt 2 ) + (-4872 - 448vr 2 ) s 2 sx + (2352 + 224tt 2 ) s 2 si + 
(-168 - 16tt 2 ) si + (120 - 16vr 2 ) s\ + (216 + 16vr 2 ) s 2 + (-168 - 16tt 2 ) s 2 

48s 2 (1 - si) (1 - s 2 ) 2 (6s 2 s? - 6s? - lls 2 s 2 + 15s 2 + 3s 2 si - 9s x + 2) x 
log(l - Si) 

24 (1 - si) (1 - s 2 ) (30s^s? - 64s 2 s? + 41s 2 s? - 7s? - 60s^s 2 + 128s 2 s 2 - 
82s 2 s 2 + 37s 2 si - 78s 2 si + 76s 2 si - 7s x - 7s\ - 7s 2 ) log (s x ) 

-48 (1 - sx) (1 - s 2 ) (s 2 s^ - s 2 sf - 5s^s? + 9s 2 s? - 5s 2 s? + s? + 9s^s 2 - 
20s 2 s 2 + 13s 2 s 2 - 5s 2 si + 12s 2 si - 12s 2 si + si + si + s 2 ) log 2 (si) 

96 (1 - sx) (1 - s 2 ) (6s^s? - 12s 2 s? + 7s 2 s? - s? - 12s^s 2 + 24s^s 2 - 
14s 2 s 2 + 7s^S! - 14s 2 S! + 12s 2 sx - s x - s 3 2 - s 2 ) log (s x ) log (s 2 ) 

-96 (1 - Si) (1 - s 2 ) (6s 3 2 sl - 12s 2 2 sl + 7s 2 s? - s? - 12s^s 2 + 24s 2 s 2 - 
14s 2 s 2 + 7S2S1 - 14s 2 si + 12s 2 si - si - s 2 - s 2 ) log (si + s 2 ) 

96 (1 - sx) (1 - 8 2 ) 2 (s 2 s? - s? + 2s 2 s? - 5s 2 s? + 3s? + s^s? - 5s 2 s?+ 
8s 2 s? - 2s? - s\s\ + 4s 2 s 2 - 4s 2 s 2 + s\ - 4s 2 si + 3s 2 si - s x - 
s 2 + s 2 ) log (1 - sx) log (si + s 2 ) 

-96 (1 - si) (1 - s 2 ) (s^s? - s 2 s? - lOs^s? + 19s^ - lls 2 s\ + 2s?- 
lls^s? + 53s^s? - 77s 2 s? + 41s 2 s? - 2s? + 21s^s 2 - 76s^s 2 + 94s 2 s 2 - 
49s 2 s 2 + 2s? - lls 2 sx + 38s^si - 46s 2 si + 25s 2 si - 2s x + s 2 - s\+ 
s 2 - s 2 ) log (sx) log (si + s 2 ) 

48 (1 - sx) (1 - s 2 ) (s 2 s? - s 2 s\ - 21s^ + 40s^sf - 22s 2 s? + 3s?- 
21s\s\ + 106s^s? - 153s 2 s? + 79s 2 s? - 3s? + s^s 2 + 40s^s 2 - 153s^s 2 + 
188s 2 s 2 - 95s 2 s 2 + 3s 2 - s^si - 22sjsx + 79s|si - 95s 2 2 sx + 50s 2 si- 
3si + 3s? - 3s 2 + 3s 2 - 3s 2 ) log 2 (si + s 2 ) 

96si (1 - 8 2 ) 2 (s 2 sf - s? + s 2 s? - 4s 2 s? + 3s? - 5s 2 s? + 8s 2 s? - 2s?+ 
7s 2 si - lls 2 si + s x - 2s\ + 5s 2 - l) Li 2 (sx) 
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/ 16 = 96(l-S 1 )(l-S 2 )(s^4_ 2s2S 4 + s 4 + 8s 3 s 3_ 17s 2 s 3 + 12s2S 3_ 

3sl + s A 2 s\ - llslsl + 32s^ - 2Qs 2 s\ - 2s^i + 12s^i - 20s^si+ 
20s 2 si — 2si + s 2 — 3s 2 — 2s 2 ) Li 2 (1 — sj — s 2 ) 



/5 = /2O1 ^ S 2 ) / 6 = / 3 (si <r* S 2 ) / 8 = / 4 (si ^ S 2 ) 

/l2 = Ao(Sl S2) /l3 = -H- S 2 ) /i 7 = /i 5 (si <B- S 2 ) 

The order a s correction dT$ / (dsids 2 ) in Eq. (15.11) to the double differential decay width 
for b — > X s 77 is the main result of our paper. 



6 Some numerical illustrations 



In the previous sections we calculated the virtual- and bremsstrahlung QCD corrections 
which were the missing ingredient in order to obtain the (Oj, O7) contribution to the double 
differential decay width for B — > X s ^ at NLL precision. The Wilson coefficient CV ie //(^) 
at the low scale (// ~ nib) which is needed up to order a s , i.e., 

CW/M = Ctz/M + ^ Cl ff (n) (6.1) 

is known for a long time (see ref. [8] and references therein). Numerical values for the input 
parameters and for this Wilson coefficient at various values for the scale /1, together with 
the numerical values of a s (/i), are given in Table [1] and Table [2j respectively. The NLL 



Parameter 


Value 


m b 


4.8 GeV 


m t 


175 GeV 


mw 


80.4 GeV 


m z 


91.19 GeV 


Gp 


1.16637 x 1(T 5 GeV~ 2 


v tb v* 


0.04 


a- 1 


137 


a s (M z ) 


0.119 



Table 1: Values of the relevant input parameters 



prediction reads 

dr 77 drg? dr# 

d^dvr : s^ 4 d^- 2 (6 - 2) 

where the first- and second term of the r.h.s. are given in eqs. (12. 5p and (15. ip . respectively. 
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jj, = mw 


0.1213 


-0.1957 


-2.3835 


fi = 2 m h 


0.1818 


-0.2796 


-0.1788 


fi = rrib 


0.2175 


-0.3142 


0.4728 


ft = m b /2 


0.2714 


-0.3556 


1.0794 



Table 2: a s (/i) and the Wilson coefficient CV ie //(/i) at different values of the scale /x 



To illustrate our results, we first rewrite the MS mass mj(/i) in eq. (16. 2p in terms of the 
pole mass m&, using the one- loop relation 



m b (fi) = m h 



^ ^8 log — + — 



Air \ rrib 3 



We then insert CV ie //(/i) in the expanded form (16 .ip and expand the resulting expression 
for dTji I \ds\ds<i) w.r.t. a s , discarding terms of order a 2 s . This defines the NLL result. The 
corresponding LL result is obtained by also discarding the order a\ term. In Fig. [6] the LL- 
and the NLL- result is shown by the short-dashed- and the solid line, respectively. 

In our procedure the NLL corrections have three sources: (a) a s corrections to the Wilson 
coefficient C 7ie //(/i), (b) expressing m fe (/i) in terms of the pole mass m b and (c) virtual- and 
real- order a s corrections to the matrix elements. To illustrate the effect of source (c), which 
is worked out for the first time in this paper, we show in Fig. [6] (by the long-dashed line) 
the (partial) NLL result in which source (c) is switched off. We conclude that the effect (c) 
is roughly of equal importance as the combined effects of (a) and (b). 

For completeness we show in this figure (by the dash-dotted line) also the result when 
QCD is completely switched off, which amounts to put /i = my/ in the LL result. 

From Fig. [6] we see that the NLL results are substantially smaller (typically by 50% or 
slightly more) than those at LL precision, which is also the case when choosing other values 
for s 2 . 
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0.2 0.4 0.6 0.8 1 



Figure 6: Double differential decay width dT^/(dsids2) as a function of s± for s 2 fixed 
at s 2 = 0.2. The dash-dotted, the short-dashed and the solid line shows the result when 
neglecting QCD-effects, the LL result and the NLL result, respectively. The long-dashed 
line represents the (partial) NLL result in which the virtual- and bremsstrahlung corrections 
worked out in this paper are switched off (see text for more details). In the frames 1), 2) and 
3) the renormalization scale is chosen to be \i — m&/2, \i = mj and fi = 2mt), respectively. 



7 Technical details about our calculations 

We first describe the general setup of our calculations and then discuss in the subsections 17.11 
and !7.2l the calculation of the virtual- and the bremsstrahlung corrections for the interference 
diagrams shown in the last lines of Fig. [3] and Fig. |5l respectively. 

The starting point is the general expression for the total decay width of the massive 
b quark with momentum pt, decaying into 3 < n < 4 massless final-state particles with 
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momenta fcj, 




where the squared Feynman amplitude \M n \ 2 is always understood to be summed over final 
spin-, polarization- and color states, and averaged over the spin directions and colors of the 
decaying b quark. It also includes a factor of 1/2 for the two identical particles in the final 
state, i.e. the photons. Furthermore, d = 4 — 2e denotes the space-time dimension that we 
use to regulate the ultraviolet, infrared and collinear singularities. 

The double differential decay rate dT 77 /(dsids2) is obtained from eq. ( 17.1 p by mul- 
tiplying the integrand on the r.h.s. with the delta functions 5 (s\ — (p& — qi) 2 /ml) and 
$ ( s 2 — (Pb — <?2) 2 / m b) [32], where pb and q\, q<i denote the four momenta of the b quark 
and the photons, respectively. For the bremsstrahlung corrections, as mentioned in sec- 
tion HI we need to consider also the triple differential decay width dT 77 / \ds1ds2ds3), where 
s 3 = (Pb ~ Qi ~ Q2) 2 / m b i s the normalized hadronic mass squared. The triple differential 
decay width is obtained by multiplying the integrand with the additional delta-function 
S (S3 — (pb — qi — q2) 2 /ml). Finally, the delta functions just mentioned and all of the delta 
functions present in eq. (17. ip can be rewritten as differences of propagators as follows [33|l3l]. 

6 - m2 )=T~- ( 2 \ n " 2 • ( 7 - 2 ) 

2m \q z — m z — i\} q z — m z + t(j J 

In this step the phase-space integrations are converted into loop integrations (which can be 
combined with possible loop integrations already present in |M n | 2 ). By subsequently doing 
tensor reductions, the (differential) decay width can be written as a linear combination 
of scalar integrals. The systematic Laporta algorithm [35], based on integration- by-part 
techniques first proposed in [36JI37], can then be applied to reduce the scalar integrals to a 
small number of simpler integrals, usually referred to as the master integrals (Mis). For this 
reduction we used the AIR and FIRE implementations [381139] of the Laporta algorithm. 
After the reduction process, it usually happens that some Mis contain propagators which 
were introduced via (17. 2p with zero or negative power. In this case the ±z0 prescription 
becomes irrelevant and as a consequence these Mis are zero. In the remaining Mis we 
convert the propagators introduced via (17. 2p back to delta-functions. Thus, we are left with 
phase-space Mis (which can contain loop integrations as well). The final task is then to 
calculate these Mis, i.e. to perform possible loop integrations together with the phase-space 
integrations. 

Very often we had to deal with Mis which we were not able to evaluate by direct in- 
tegration of their integral representation in terms of Feynman parameters and/or phase- 
space parameters. A powerful tool to be used in these cases is the differential equation 
method [33l [3^1 SD1 Si] - The goal of this method is to employ the output of the reduction 
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procedure for a given topology to build differential equations which are satisfied by the Mis 
of that topology. In our case, we consider differential equations w.r.t. sx and s 2 and also 
w.r.t. S3 for case of bremsstrahlung corrections. With these methods we were able to ob- 
tain analytic expressions for all master integrals appearing in the calculation of the various 
diagrams. 



7.1 Details about the calculation of virtual corrections 

To illustrate our methods for the virtual corrections, we take as an example the interference 
diagram shown on the last line of Figure [3j In this case we have five Mis. Four of them can 
be solved by means of direct integration on Feynman parameters. To calculate the last one 
(called Pun), we solve the differential equations w.r.t. si and s 2 and get the solution which 
we denote as Qnn- At this level Qun contains integration constants (which are not fixed 
by the differential equations). To get the integration constants, we proceed in the following 
way. Using Feynman parametrization for the loop integral, we write the MI P mi as 

Pun = Si e S2 e (l - si - s 2 )~ e / g (si,s 2 ,e,u,v,y) dudvdy , (7.3) 

Jo 

where u, v and y are Feynman parameters (all of them running from to 1). The factor 
S V S 2 € (^ ~ s i ~ s 2)~ t is coming from the phase-space (see eq. HB.2[) in Appendix IB.ip . One 
then can put s 2 = in go(si, s 2 , e, u, v, y) and integrate on the Feynman parameters, which 
defines the new function 

9i{si,e) = g (s 1 ,0,e,u,v,y)dudvdy . (7.4) 
Jo 

We managed to work out the leading term of the expansion of <7i(si,e) on S\ around zero, 
which is proportional to sj~ 2 . From this, we immediately get the leading term of the expansion 
of Pun on s 2 and si (which is proportional to s^/sf). Comparing the result of this calculation 
with the corresponding expansion of the solution Qnn of the differential equations, we could 
determine all integration constants. 



7.2 Details about the calculation of bremsstrahlung corrections 

To illustrate our methods for the bremsstrahlung corrections, we take as an example the 
interference diagram shown on the last line of Figure |5j In this case we obtain three Mis, 
denoted by Poo, P10 an d Pn- Writing the diagram as a linear combination of the Mis, we 
see that the leading power (w.r.t. S3) of all three coefficients (in front of the Mis) is of 
the order l/s 3 . Keeping in mind that we are taking into account only terms proportional 
to the leading power in S3 at the level of the triple differential decay width (as extensively 
discussed in section HJ), it is sufficient to work out the Mis to zero-th power in S3, including 
the epsilontic regulator, i.e., s®~ ne (in our case only n = 1 and n = 2 occur). 

The simplest MI, P o, which corresponds to the pure phase-space (see eq. (IB. 4ft in 



Appendix IB.2p . can be easily solved by means of direct integration. For P 10 the solution of 
the differential equation w.r.t. s 3 can be represented in the limit s 3 — > in the form 

P10 = s 2 , e) 53"' + a 2 (si, s 2 , e) s^ 2e , (7.5) 
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where the function ai(si, s 2 , e) is fully determined. To find the function a 2 (si, s 2 , e), we use 
differential equations w.r.t. si and s 2 . In this way, we could find a 2 (si, s 2 , e) up to integration 
constants. To determine these constants, we managed to calculate the MI for specific values 
of si, s 2 and S3 — > 0. In the same way we also calculated the MI P n . 



8 Summary 

In the present work we calculated the set of the 0(a s ) corrections to the decay process 
B — > X s 77 originating from diagrams involving the electromagnetic dipole operator O7. To 
perform this calculation it is necessary to work out diagrams with three particles (s-quark 
and two photons) and four particles (s-quark, two photons and a gluon) in the final state. 
From the technical point of view, the calculation was made possible by the use of the Laporta 
Algorithm to identify the needed Master Integrals and by applying the differential equation 
method to solve the Master Integrals. When calculating the bremsstrahlung corrections, we 
take into account only terms proportional to the leading power of the hadronic mass. We 
find that the infrared and collinear singularities cancel when combining the above mentioned 
approximated version of bremsstrahlung corrections with the virtual corrections. The nu- 
merical impact of the NLL corrections is not small: for dT^/(dsi g?s 2 ) the NLL results are 
approximately 50% smaller than those at LL precision. 
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A Explicit results for the functions V; L defining the vir- 
tual corrections 

The functions V{ appearing in eq. (j3.5p read 

vt = (1-s 1 -s 2 )(-16tv 2 sIs 5 1 -16it 2 sIsI + A8it 2 s 2 2 s 5 1 +4:87i 2 sIs 2 1 -4:87i 2 S2sI- 

48vr 2 s^i + 16vr 2 ^ + 16n 2 s 5 2 + (2112 - 104vr 2 ) s\s{ + (392tt 2 - 6384) s 3 s\+ 

(6672 - 532vr 2 ) s 2 2 s\ + (288vr 2 - 2736) s 2 s\ + (336 - 60vr 2 ) s\+ 

(392tt 2 - 6384) s A 2 s\ + (19584 - 1224vr 2 ) s\s\ + (l452vr 2 - 20784) s 2 2 s\+ 

(7872 - 600tt 2 ) s 2 s\ + (44vr 2 - 288) s\ + (6672 - 532tt 2 ) s^s 2 + 

(1452tt 2 - 20784) s 3 s\ + (23904 - 1728vr 2 ) s 2 2 s\ + (740tt 2 - 10128) s 2 sj+ 

(336 - 28tt 2 ) si + (288tt 2 - 2736) sfsi + (7872 - 600tt 2 ) s 3 2Sl + 

(740tt 2 - 10128) s 2 Sl + (5376 - 392tt 2 ) s 2Sl + (28tt 2 - 384) 8l + 

(336 - 60tt 2 ) s\ + (44tt 2 - 288) s\ + (336 - 28vr 2 ) s 2 2 + (28tt 2 - 384) s 2 ) 

V2 = _96(l -si)(l-s 2 )(l-si -s 2 ) (3s^s? - 4s^s? + s 2 s? - 5s^s 2 + 7s 2 s 2 - 
2s 2 s 2 — s 2 + 2S2S1 — 3s 2 Si + 3s 2 s x — s 2 ) log (si) 

v 3 = -24 (1 - si) (1 - s 2 ) (1 - si - s 2 ) (2s 2 s? - 2s 2 s\ - As 3 s 3 + 6s 2 2 s 3 - 
3s 2 sl + sl + 6s 3 2 sl - 16s 2 s 2 + 12s 2 s 2 - 3s|si + lOs 2 ,^ - 12s 2 si+ 
81 + s 3 2 + s 2 ) log 2 ( Sl ) 

w 4 = 48 (1 - si) (1 - s 2 ) (1 - si - s 2 ) (6s^s? - 12s 2 ^ + 7s 2 s? - s? - 12s|s?+ 
24s 2 s 2 - 14s 2 s 2 + 7S2S1 - 14s 2 si + 12s 2 si - si - s| - s 2 ) log (si) log (s 2 ) 

w 5 = 48 (1 - si) (1 - s 2 ) (1 - si - s 2 ) (6sijs? - 12s 2 s? + 7s 2 s 3 - s 3 - I2s 3 sl+ 
24s 2 s 2 - 14s 2 s 2 + 7sls 1 - 14s 2 ,si + 12s 2 Si - s x - s 2 - s 2 ) log (si) x 
log(l - si - s 2 ) 

w 6 = -96(l-si) 2 (l-s 2 ) 2 s 2 (st + 2s 2 s?-2s? + s 2 s 2 -4s 2 s 2 + s 2 -2s 2 Sl + 
3s 2 si - 2si + si + l) log (si) log (si + s 2 ) 

w 7 = 48(l-s 1 )(s 2 -l) 2 s 2 (l-s 1 -s 2 ) (6s 2 s?-6s?-lls 2 s 2 + 15s 2 + 
3s 2 si-9si + 2)log(l-si) 

w 8 = 96(l-si)(s 2 -l) 2 (s 2 s?-s? + 2s 2 4-5s 2 4 + 3st + s^?-5s 2 s3+ 
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8s 2 s? - 2sf - s^s 2 + 4s 2 s 2 - As 2 s\ + s\ - As 2 Si + 3s 2 s 1 - s 1 - s% + s 2 ) x 
log(l - si)log(si + s 2 ) 

v 9 = 48 (1 - si) (1 - s 2 ) (sjsl - s 2 sl - 9s 3 2 sf + Wsjs* - 8s 2 s? + s{ - 9s^?+ 
46s 3 2 sl - 67sjsl + 35s 2 sl -s\ + s\s\ + 16s^s? - 67s 3 2 sl + 84s 2 ,s 2 - 
43s 2 s 2 + s 2 — s^Si — 8S2S1 + 35S2S1 — 43S2S1 + 22s 2 Si — si + s\ — s\+ 
si - s 2 ) log 2 (si + s 2 ) 

U10 = -96 (1 - si) (1 - s 2 ) (1 - si - s 2 ) (s^s? - s 2 s? + s 2 s 2 - 3s 2 s 2 + 2s 2 s 2 - 
s 2 - sl-sx + 2^2^! + S 2 S X - sf) log (1 - Si - s 2 ) 

v u = 24 (1 - si) (1 - 8 2 ) (l-st- 8 2 ) (Gs 3 2 sl - 12s 2 s? + 7s 2 s? - sf - 12s 3 2 sl+ 
24s 2 s 2 - 14s 2 s 2 + Ts^i - 14s 2 si + 12s 2 si - si - s 2 - s 2 ) x 
log 2 (1 - si - s 2 ) 

Vl2 = 96 Sl (l- S 2) 2 (l-Si-s 2 )(s 2 s^-sf + s^?-4s 2 s? + 3s?-5s^ 2 + 
8s 2 s 2 — 2s 2 + 7s 2 si — lls 2 si + si — 2s 2 + 5s 2 — l) Li 2 (s x ) 

via = 96(l-si)(l-s 2 )(l-si-s 2 ) (s 2 s? - 2s 2 s* + s* + 85^ - 17s 2 s?+ 
12s 2 s? - 3s? + s\sl - I7s 3 sl + 32s 2 s 2 - 20s 2 s 2 - 2s|si + 12s^ Sl - 
2OS2S1 + 20s 2 si - 2si + s 2 - 3s2 - 2s 2 ) Li 2 (1 - s x - s 2 ) 

Vl4 = V2 ( Sl o s 2 ) V 15 = V 3 (St «-» s 2 ) v w = u 5 (si «-»> s 2 ) 

vn = v 6 (s 1 <^s 2 ) v 18 = v 7 (si <H> s 2 ) v 19 = u 8 (si O s 2 ) Vjjo = Vufai s 2 ) 



B Relevant phase- space formulas 

B.l Double differential phase-space for the 3-particle final state 

The kinematical variables Si and s 2 are defined as 

81 = 2 — ; s 2 = 2 — > R1 

where p& and ^ denote the four-momenta of the 6-quark and the photons, respectively. The 
kinematics of the process b — > S77 is fully described by Si and s 2 . The formula for double 
differential decay width is therefore free of additional phase-space integration variables. It 
reads 

The variables Si and s 2 vary in the range 
< si < 1 ; < s 2 < 1 - si . 
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B.2 Triple differential phase-space for the 4-particle final state 

The kinematical variables si, s 2 and S3 are defined as 

(Pb ~ qi? {Pb - g 2 ) 2 (Ps + kf 

Si — g j s 2 — 2 ' S 3 — o \ a - 6 ) 

mi m^ mi 

where pi,, p s , k and g« denote the four-momenta of the 6-quark, the s-quark, the gluon and 
the photons, respectively. The kinematics is fully described in terms of five phase-space 
variables x\, x 2 , ^3, ^4 and £5 as given explicitly in eqs. (3.6), (3.9) and (3.10) in ref. [12]. 
By identifying hi, k2, k% and k$ given there with the four- momenta of the two photons, the 
s-quark and the gluon, respectively, we easily derive from the information in [12] the formula 
for the triple differential decay width. We remind the reader that in this paper we consider 
only the range in s\ and S2 with 

< si < 1 ; < s 2 < 1 - si , 
which is also accessible to the three-body decay b — > S77. For this case we obtain 

™ b s 3 ( s i*2 - s 3 y e (i- Sl -s 2 + s 3 y e x 



ds 1 ds 2 dsz (1 - 2e)T 2 [l - 2e] 

J dx, dx 5 [x 4 (1 - x 4 )]" e [x 5 (1 - x 5 )]" 1/2 - £ |M 4 | 2 , (B.4) 

where x±, x 2 and x 3 (appearing in |M 4 | 2 ) are understood to be expressed in terms of si, s 2 
and s 3 according to 

s 3 sis 2 - s 3 
Xi = si ; x 2 = — ; x 3 = r . (B.5) 

Si (1 - Si) (S! - S3) 

x 4 and x 5 vary between zero and one, while s 3 G [0, Si s 2 ]. 



C Renormalization constants 

In this appendix, we collect the explicit expressions of the renormalization constants needed 
for the ultraviolet renormalization in our calculation (see section [3]). 

The operator (9 7 , as well as the 6-quark mass contained in this operator are renormalized in 
the MS scheme 



zW=l + ^^ + 0(a>) ; ^=l-^^ + 0(« 2 ). (C.l) 

All the remaining fields and parameters are renormalized in the on-shell scheme. The on-shell 
renormalization constant for the 6-quark mass is given by 
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while the renormalization constants for the s- and fc-quark fields are 



3 — 2e / /i \ 2e a s (/x) 



The various quantities appearing in section [3] are defined to be SZ = Z — 1. 
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